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On Mars, atomic oxygen controls the carbon dioxide radiative cooling of 


the upper atmosphere and the presence of an ozone layer near the poles. 

To remotely probe meridional transport of O atoms from the summer to the 
winter hemisphere and the descending flow in the winter polar regions, the 
O, Herzberg Il atmospheric emission could be used as a proxy. This emission 
is quite weak on Earth’s nightside, but itis prominent in the Venus night 
airglow, and it has not previously been observed on Mars. Here we report the 
limb detection of the O, Herzberg II visible bands in the Mars nightglow with 
the NOMAD ultraviolet—visible spectrometer onboard the European Space 
Agency’s Trace Gas Orbiter. The emission layer reaches up to hundreds of 
kilorayleighs in the limb viewing geometry. It is mainly located between 

40 kmand 60 kmat high latitudes during the winter season, consistent with 
three-body recombination of oxygen atoms. This O, nightglow should be 
observable from a Martian orbiter as well as from the Martian surface with 
the naked eye under clear sky conditions. These observations pave the way 
to future global observations of the Martian atmospheric circulation with 
simpler lower-cost instrumentation. 


The mesosphere of Venus and Mars is still poorly known as data are 
scarce and global circulation models can hardly reproduce them. Its 
composition, circulation pattern and variability are still not well under- 
stood. It is important to characterize the Martian upper atmosphere 
to determine the structure of the terrestrial planets’ atmospheres and 
understand their different evolution. Observations of airglow emissions 
is a powerful remote-sensing tool that provides insight into the density 
variations and global transport in this region. 

The Mars airglow emissions on the dayside have been extensively 
studied based on observations performed with instruments onboard 
the Mariner, Mars Express (MEx), Mars Atmosphere and Volatile Evo- 
lution Mission (MAVEN), Trace Gas Orbiter (TGO) and Emirates Mars 
Mission (EMM). Following the first observations from the Mariner 6 and 


7 flybys’, later studies have described the distribution of the ultraviolet 
CO Cameron, CO,* UV doublet and Fox-Duffendack-Barker bands, 
HLyman-a, 01130 nm and 136 nm multiplets, and N, Vegard-Kaplan 
emissions, and their sources and dependence on latitude, solar lon- 
gitude and solar activity. Recently, observations of the [O 1] 557.7 nm 
and the 630-636 nm forbidden emissions made with the TGO Nadir 
and Occultation for Mars Discovery (NOMAD) instrument and their 
seasonal variability have been reported’. 

On the nightside, only the ultraviolet and infrared emissions have 
been explored with spectral instruments onboard the MEx and MAVEN 
missions. The middle ultraviolet spectrum is dominated by the v’ = 06 
and y bands of nitric oxide excited by radiative association of nitrogen 
and oxygen atoms. Although this emission is present at all latitudes 
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and local times, extensive mapping has shown that it is enhanced at 
high winter latitudes in both hemispheres’. This brightening is a con- 
sequence of the transport of O and N atoms from the summer to the 
winter hemisphere where they are carried downwards by vertical winds 
and diffusion and recombine in the 40-60 km region. The O, infrared 
atmospheric band at 1.27 pm’ ’ and the OH Meinel emission at 1.45 um 
and 2.9 um” have also been measured in the Martian polar winter 
atmosphere over a similar range of altitudes. They also result from 
meridional transport and winter polar descent of O atoms produced 
by photodissociation of CO, in the summer hemisphere. 

No nightglow spectrum has been observed so far in the visible 
domain. Scattered solar light hinders such measurements from Earth’s 
ground and no space-borne instrument has observed the Mars visible 
nightglow. In this Article, we describe new observations made with the 
ultraviolet and visible (UVIS) channel, a spectrometer of the NOMAD 
instrument flying on the ExoMars TGO. They cover both hemispheres, 
but their number is limited by operational constraints. Nevertheless, 
visible nightglow emissions have been observed on several TGO orbits. 
The observations reported here were episodically collected between 
March 2020 and October 2022 with the UVIS instrument. The TGO 
satellite has been orbiting Mars since October 2016 ona circular orbit 
at ~400 km, inclined 74° relative to the equator with an orbital period 
of -2 h. The UVIS instrument covers the 200-650 nm spectral range 
and can observe in both the nadir and the solar occultation modes. 
The UVIS spectral resolution varies from 1.2 nm at 200 nm to 1.6 nm 
at 650 nm™”. For the observations described in this study, the TGO 
spacecraft was oriented so that the line of sight of the nadir channel was 
pointing in the limb direction”. Two pointing modes were used during 
this period. Inthe first one, ‘inertial pointing’, the line of sight was kept 
fixed in the inertial space so that the line of sight scanned a range of 
altitudes usually from 300 km, down to 40 km, then up to 300 km. In 
the second one, most frequently used in this study, the areoid altitude 
of the tangent point in the atmosphere remained within ~20 km of a 
pre-set value while different latitudes were probed due to the motion 
of the spacecraft. The spectra were integrated over 20 s. The resolution 
at the tangent point, taking into account the UVIS field of view during 
the spectral integration provided a vertical resolution of 10 +5 km. 

The NOMAD-UVIS instrument occasionally collected limb spectra 
during nightside conditions, with solar zenith angles greater than 
105° at the tangent point of the line of sight. A band structure clearly 
appeared in some individual bright spectra between 400 nm and 
600 nm. The wavelengths and relative intensities indicate coincidence 
with the strongest bands of the v’ = 0 progression of the O,c'Z, > X°x, 
Herzberg Il transition”. 

A total of 1,228 spectra were acquired between 20 km and 
70 km and with a solar zenith angle >105°, out of which 122 were kept 
(Methods). The average spectrum is shown in Fig. 1. The (0,5) to (0,11) 
Herzberg II bands are clearly identified with relative intensities similar 
to those observed on Venus” and in the laboratory. Beyond 620 nm, 
the noise level increases, making detection of the (0,12) band uncertain. 
The total intensity for bands v’ = 5 to v’ = 12 derived from the spec- 
tral fit combined with the relative band intensity in the laboratory 
is 108 + 11 kR. The estimated error combines contributions from the 
absolute calibration and measurement noise including the uncertainty 
in the spectral baseline. 

The confirmed spectral detections spread over 13 orbits are shown 
in Fig. 2. A first group of detections extends from 50° S to 81° S when 
the solar longitude L, varied from 94° to 158°. A second one covers 
the L, range 225-330° between 65° N and 81° N. The mean spectrum 
is shown in Fig. 1. 

The observed wavelengths perfectly match those in the laboratory 
and Venus Herzberg II spectra. The small differences in relative band 
brightness in Fig. 1 result from the presence of the background noise 
level in the UVIS data. The relative band intensities in the brightest 
individual spectra show excellent agreement of the spectral shape. 
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Fig. 1| Average high-latitude nightglow spectrum. The identifications of the 
v’ = O Herzberg II bands are indicated by the vertical dashed lines. The purple 
curve is the reference Venus nightglow spectrum scaled to the UVIS mean 
spectrum. A total of 122 high-latitude NOMAD-UVIS spectra have been averaged 
(black curve). 


There is no indication of a departure of the intensity distribution com- 
pared with the Venus or laboratory spectra. By analogy with the Venus 
nightglow spectrum, the presence of other weak O, band systems, such 
as the Chamberlain A?A, > a'A,, Herzberg 1 A*Z, > X°Z, andc'Z, > a'A, 
bands, may be expected in the region of the (0-6) to (0-12) Herzberg 
Il bands. On Venus, their contribution is about 7%, 5% and 4% of the 
Herzberg Il intensity, respectively”. However, the signal-to-noise ratio 
of the UVIS instrument is clearly too low to confirm their presence in 
the UVIS nightglow spectra. 

In Fig. 2, inertial limb detections are represented with dots and 
limb tracking observations are shown with triangles. Each dot colour 
corresponds to one TGO orbit. Grey symbols represent observations 
without detection. Figure 2a illustrates the concentration of the detec- 
tions under high-latitude winter conditions. Figure 2b indicates that 
the visible nightglowis observed between 34 km and 66 km, generally 
polewards of 60°. The statistics of the detections are too small to obtain 
a statistically significant mean limb profile. 

Figure 3a shows the histogram of the altitudes of the spectra 
binned into 1 km intervals corresponding to a positive detection at 
the linear Pearson correlation coefficient R > 0.32 level. Amaximum 
of 35 detections is observed at 42 km. The brightest single spectrum 
reached an intensity of 290 kR at 41 km and 76° N. The detection lower 
limit for single 15 s spectra is on the order of 20 + 10 kR. Figure 3b pre- 
sents the histogram of the detection relative frequency of occurrence 
versus altitude. They are mainly grouped between 35 km and 59 km. The 
highest value is also observed at 42 km where the occurrence is a large 
as 76%. The detection frequencies at higher altitudes remain below10%. 

The detection of the visible Herzberg II bands in the Mars dayglow 
makes its measurement a promising technique to derive the atomic 
oxygen density and its variations in the middle atmosphere. It is pro- 
duced by three-body recombination of atomic oxygen, following the 
scheme: 


0+0+M—>05+M (1) 


where M is the atmospheric number density. This is the accepted 
mechanism for the excitation of the O, nightglow in atmospheres”. 
Models based solely on three-body recombination of O atoms have 
been shown to adequately explain the observed O, Herzberg Il intensity 
in the Venus nightglow’ and the O,'A nightglow on Mars. Other pro- 
cesses have been mentioned as potential sources of O, c'Z molecules 
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Fig. 2 | Distribution of detections of the O, Herzberg II nightglow. 
a, Latitude and season of the detections. b, Areoid altitude-latitude map of the 
observations. The triangles indicate limb tracking and the dots indicate inertial 
pointing orbits. Each dot colour corresponds to one TGO orbit. Small grey dots 
correspond to observations without a measurable Herzberg II spectral signature. 


on Venus, but they have been considered negligible in models of the 
O Venus and Mars nightglow'*”. In the absence of evidence for other 
source processes, we consider process (1) as the only important source 
of O,c'Z excitation. 

Model simulations for a zonally averaged atmosphere have been 
performed for different L, and latitudes in both hemispheres. The vol- 
ume emission rate is calculated based on process (1) and the relevant 
Mars Climate Database (MCD) temperature and density distributions. 
Figure 4 shows the calculated limb profile of the brightness of the 
Herzberg II bands. Slightly lower values are predicted at latitudes 
beyond 75° but they are expected to drop considerably at lower lati- 
tudes. The calculated peak brightness decreases to a few kilorayleighs 
at 60° N and moves up to 58 km. For the sake of completeness, the 
calculated limb profile for 60° N is also shown, although no night- 
glow observation has been performed in the limited database. Some 
intensity measurements from individual spectra are also shown for 
comparison. They either closely agree with the model value (79° N, 
76° N, 76° S) or exceed it by a factor of 2 (75° N). The measurement at 
L,=50°, 65° Shas been collected in May 2023 and shows a discrepancy 
relative to the model by about an order of magnitude. This indicates 
that the O density was significantly higher than the value predicted by 
the MCD at the time of the observations. This difference is observed at 
a moderately high latitude, suggesting that the modelled meridional 
transport does not provide enough O atoms at these lower latitudes. 
If discrepancies are demonstrated to be reliable, they should provide 
new insights into the dynamical processes inthe atmosphere, thereby 
improving the reliability of three-dimensional GCMs. Additional obser- 
vations are necessary to confirm whether the discrepancies at these 
latitudes are common or occasional. 

The clear spectral signature and its analogy with the laboratory 
and Venus nightglow, together with the general agreement of the 
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Fig. 3 | Altitude distribution of the nightside limb observations. 
a, NOMAD-UVIS observations (light blue) and positive detections (dark blue). 
b, Rate of occurrence of detections. 


predicted altitude and brightness, confirm that the observed spectrum 
is produced by three-body recombination of O atoms. The enhanced 
brightness at high latitudes is the signature of the global transport 
ascending from the sunlit summer polar regions and descending to 
the dark winter hemisphere. 

The O, a'4A, > X°Z, infrared atmospheric emission at 1.27 yum is 
another recombination emission also observed in the Martian night- 
glow. The peak altitude of the emitting layer varies between 40 km 
and 65 km reaching up to 12 MR at 40 km in both winter regions at 
latitudes beyond -70°%°. In the case of the a'A, metastable state, the 
recombination efficiency of the O atoms is on the order of 70% and 
collisional deactivation is negligible at the altitude of the emitting 
layer. Consequently, the intensity of the infrared atmospheric system 
is about 60 times higher than the Herzberg II bands. 

Collisions between O and CO, in the thermosphere populate the 
vibrational levels of CO,, resulting in enhanced radiation of the 15 um 
emission to space and atmospheric cooling”. Consequently, the O 
density governs the efficiency of the non-local thermodynamic equi- 
librium (non-LTE) radiative cooling in the upper atmosphere. Although 
this process essentially occurs at higher altitude than most of the O, 
airglow emission, the O peak density derived from the nightglow verti- 
cal profile may be used to constrain the O density at higher altitudes. 
Therefore, the determination of its distribution and variability is of 
prime importance in understanding and modelling the thermal struc- 
ture of the Martian upper atmosphere. Transport of O atoms across 
the terminator also has a key role asa precursor inthe formation of an 
ozone layer near the south winter pole”. 

Future measurements will allow better determination of the emis- 
sion peak altitude, its brightness and variability. They will provide 
unique constraints on three-dimensional models of the global circu- 
lation” and dynamical processes at play in the middle atmosphere. 
By accumulating nightside limb observations, it will be possible to 
describe the mean O density distribution at high winter latitudes and 
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Fig. 4 | Model simulations of the limb brightness distribution of the O, 
Herzberg II bands. The simulated limb profiles are calculated from expression 
(2) (Methods) based on densities and temperature from the MCD for 

different latitudes and seasons. The circles indicate UVIS observations using 
corresponding colours. The observed intensity values are: 117 kR, 45 km (orange); 
117 KR, 51 km (red); 162 kR, 40 km (blue); 290 KR, 41 km (cyan); 86 kR, 54 km 
(black). The error bars are estimated at 10% of the brightness. No detection at 
latitudes as lowas 60° Nis available in the current database. 


its variability. In particular, additional observations will determine how 
accurately the transition of the bright nightglow from one hemisphere 
tothe other is predicted by three-dimensional models. The question of 
a possible hemispheric asymmetry in the spatio-temporal airglow dis- 
tribution will also be examined based ona larger number of detections. 

The brightness of the O, Herzberg bands of several hundred 
kilorayleighs in the visible opens the possibility to observe them in 
the Martian high-latitude regions with simple space-borne cameras 
in the visible domain. We also note that the sensitivity threshold of 
the human eye in the green for nightside vision is estimated on the 
order of 1 kR”, decreasing at shorter and longer wavelengths”. For the 
Herzberg bands, which are widespread over most of the visible range, 
a likely wavelength-averaged threshold value is on the order of a few 
kilorayleighs, much lower than the bright emission described in this 
report. The polar winter visible nightglow is therefore expected to be 
observable with the naked eye by future astronauts from orbit or the 
Martian ground. Incidentally, we note that the oxygen green line at 
557.7 nm is not observed in UVIS bright individual spectra nor after 
co-addition of nightglow spectra (Fig. 1). This non-detection sets an 
upper limit of about 1 kR on the line brightness. It is in full agreement 
with its absence in the laboratory recombination spectra and the lack 
of detection by the Mars 5 visible spectrometer setting an upper limit 
of 50 R” onthe nightglow limb brightness of the green line. 


Methods 

Data processing 

During the pipeline processing of the spectra, the instrumental back- 
ground and the dark current on the charge-coupled device (CCD) 
frame are first removed. In a second step, the noise is corrected from 
cosmic rays, anomalous and hot pixels. Next, the 81 fully illuminated 
CCD lines are binned to increase the signal-to-noise ratio of the indi- 
vidual spectra. Ina third step, the count rate is converted into physical 
units (kilorayleighs, 1 kR corresponds to a brightness of 10° photons 
per square centimetre per second emitted in 4m steradians), using 
laboratory measurements and in flight calibration”. In addition to the 
statistical error onthe count rate and removal of the background noise, 
a possible systematic error of -10% is associated with the uncertainties 
of the absolute instrumental calibration. 


For verification of positive spectral detections, we calculate the 
correlation coefficient between the individual NOMAD spectra and the 
Venus nightglow average spectrum”’, between 410 nm and 560 nm at 
aresolution of 3 nm. This Venus reference spectrum was obtained by 
summing nightglow spectra collected between 92 kmand100 km with 
the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS-M) 
on board Venus Express during a total of 120 h*. The Venus reference 
spectrum has been scaled to the UVIS averaged spectrum, follow- 
ing a least-squares fit. NOMAD-UVIS spectra with a linear correla- 
tion coefficient R > 0.32 have been retained and those with a too low 
signal-to-noise ratio have been removed. Application of the statistical 
Fisher test indicates that this correlation coefficient is significantly 
different from zero with a confidence level of 99%. A subsequent visual 
examination identified a number of contaminated spectra that have 
been removed from the database. 


Model simulations 
The volume emission rate of the band system based on process (1) is 
given by: 

A 


T 2 
=T M Ta] 2) 


where [M] is the total number density, and [O] and [CO,] are the O and 
CO, densities from the MCD; € = 0.03 is the effective yield of the c'Z, 
state including cascades from higher lying states; A = 0.17 s tis the total 
transition probability of the c-X transition, that is, from (c, v’ = 0) toall 
(X, v’); K=7.5 x 10” (7/300) > is the total recombination coefficient 
of the three-body recombination of O in CO,; keg) = 3.1 x 10™ cm? st 
and ko = 5.9 x 10” cm“ s” are the quenching coefficients of the c£, 
state by CO, and O, respectively. 

We adopt values deduced from the studies of the Venus airglow” 
based on laboratory measurements” and Venus Express observa- 
tions”°’’. They have been slightly adapted to match the O, nightglow 
limb observations with the VIRTIS spectral imager**”’. The K coefficient 
and its temperature dependence are adopted from ref. 30. The recom- 
bination rate is multiplied by 2.5 to account for the higher efficiency of 
CO, relative to N, as a third body”. For comparison with our observa- 
tions, the O and CO, densities are taken from the MCD version 6.1” for 
minimum solar activity conditions. Limb profiles are calculated by 
integrating the volume emission rate along the line of sight. 


Data availability 

The NOMAD-UVIS spectra can be downloaded from ESA’s SA 
archives at https://archives.esac.esa.int/psa/#! Table%20View/ 
NOMAD=instrument (select UVIS from the list of instruments and 
‘Level 3 Calibrated’ from the processing level). Observed limb inten- 
sities and model calculations supporting Fig. 4 are available from 
BIRA-IASB repository at https://doi.org/10.18758/71021084 or from 
the corresponding author upon reasonable request. 
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